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Jitter radiation images, spectra, and light curves from a 
relativistic spherical blastwave 
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ABSTRACT 

We consider radiation emitted by the jitter mechanism in a Blandford-McKee self- 
similar blastwave. We assume the magnetic field configuration throughout the whole 
blastwave meets the condition for the emission of jitter radiation and we compute the 
ensuing images, light curves and spectra. The calculations are performed for both a 
uniform and a wind environment. We compare our jitter results to synchrotron results. 
We show that jitter radiation produces slightly different spectra than synchrotron, in 
particular between the self-absorption and the peak frequency, where the jitter spec- 
trum is flat, while the synchrotron spectrum grows as v^^^. The spectral difference 
is reflected in the early decay slope of the light curves. We conclude that jitter and 
synchrotron afterglows can be distinguished from each other with good quality obser- 
vations. However, it is unlikely that the difference can explain the peculiar behavior 
of several recent observations, such as flat X-ray slopes and uncorrelated optical and 
X-ray behavior. 
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1 INTRODUCTION 

The synchrotron external shock model (Meszaros & Rees 
1997) has been used with great success to model GRB after- 
glows since their discovery (Wijers, Rees & Meszaros 1997; 
Piran 1999; Panaitescu & Kumar 2001; Zhang & Meszaros 
2004). However, its physical foundations have been put un- 
der debate, principally because a coherent equipartition 
magnetic field is assumed to be generated at the shock front. 
In addition, the fraction of internal energy that is given to 
the magnetic field is assumed independent of the location 
in the blastwave (the distance from the shock) and of the 
properties of the shock (e.g. its Lorentz factor). 

The two latter assumptions have been tested phe- 
nomenologically on afterglow data. Rossi & Rees (2003) con- 
sidered a magnetic field that is allowed to decay behind the 
shock front. They concluded that afterglow observations re- 
quire a magnetic field that has an approximately uniform 
intensity throughout the blastwave (or at least within a dis- 
tance R/r^ from the shock front). Yost et al. (2003) studied 
the possibility that the fraction of the internal energy stored 
in the magnetic field component grows with the decrease of 
the Lorentz factor of the fireball. They concluded that this 
effect could explain a sample of radio afterglows that have 
flat late time radio decay. 

Even though these works present interesting constraints 



on the behavior of the magnetic field, they do not address 
the mechanism that generates it, its structure, and whether 
the radiation mechanism is really synchrotron or not. An- 
alytic (Medvedev & Loeb 1999) and numerical work (Silva 
et al. 2003; Nishikawa et al. 2003; Frederiksen et al. 2004; 
Medvedev et al. 2005; Spitkovsky 2008; Chang, Spitkovsky 
& Arons 2008) have started to unveil recently the mechanism 
by which a magnetic field component is generated behind a 
relativistic coUisionless shock front. All the studies find that 
the magnetic field is generated through some kind of two- 
stream instability, such as the Weibel instability (Weibel 
1959). Magnetic fields are generated behind the shock with 
an energy of the order of 30 per cent of the total available in- 
ternal energy. They are observed to subsequently decay and, 
possibly, to stabilize at a level of a fraction of a per cent of 
the internal energy. Whether they decay even further deep 
into the shocked material is a matter of debate. The struc- 
ture of the field is observed to be highly tangled and planar 
behind the shock front and to grow more coherent deep into 
the blastwave. 

The correlation length of the magnetic field right behind 
the shock is so short that the synchrotron approximations 
do not hold and a new radiative regime sets in. This form of 
radiation, called jitter radiation (Medvedev & Loeb, 1999; 
Medvedev 2000, 2006 or, sometimes, diffusive synchrotron 
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radiation, Fleishman 2006ab), has been studied analytically, 
semi-analytically and numerically (Medvedev & Loeb 1999; 
Medvedev 2000; Frederiksen et al. 2004; Medvedev 2006; 
Fleishman 2006a; Medvedev et al. 2007; Workman et al. 
2008). In this work we apply the theory of jitter radiation 
to compute the emission from a Blandford & McKee (1976) 
blastwave, taking into account its structure, the special rel- 
ativistic effects, and the light propagation time. The results 
include spectra, lightcurves and images, analogously to what 
is done by Granot, Piran & Sari (1999a) for the synchrotron 
case. 

This paper is organized as follows: In § 2 we describe 
the numerical methods applied to integrate the radiation 
spectrum over the BM solution and over surfaces of equal 
arrival time; in § 3 we describe our results and in § 4 we 
summarize and discuss them. 



2 MODELING THE AFTERGLOW 

In this section we describe the process by which we solve 
for the physical shape, evolution, and internal properties of 
a Gamma-Ray Burst (GRB) blast wave. We loosely follow 
the outline used by Granot et al. (1999a) but we expand on 
their work by allowing for both a constant density ISM and 
a wind environment such as would possibly be found near 
the progenitor stars of long duration Gamma Ray Bursts. 
We will describe the physical setup of the problem and 
then describe the method used to compute images, spectra, 
and light curves assuming that the whole blastwave radiates 
jitter radiation. See Medvedev & Loeb (1999), Medvedev 
(2000) and Workman et al. (2008) for a thorough discussion 
of the conditions under which jitter radiation is produced. 



2.1 The Blandford McKee Solution 

To calculate the spectrum of a GRB afterglow emitting jitter 
radiation we first have to compute the local physical quanti- 
ties in the blastwave and the time evolution of the afterglow 
in the relevant geometry of the external medium. As soon as 
the external shock has swept up a mass comparable to the 
rest mass of the fireball over its Lorentz factor, the blastwave 
enters the self-similar stage (Blandford & McKee 1976). 

Following Blandford & McKee (1976) we define a simi- 
larity variable as: 



xW = l + 4(4-fc)7/(l-;^) 



(1) 



where R is the radius of the shock, r is the distance to the 
center of the shock, 7/ is the Lorentz factor of the mate- 
rial immediately behind the shock front and k describes the 
properties of the external medium. For a constant density 
profile, fc = 0, while for a wind environment fc = 2. Us- 
ing this similarity variable, we can solve for the comoving 
number density {n), the comoving energy density (e'), and 
the Lorentz factor (7) of the material anywhere within the 
expanding volume with the following equations: 



, 1 (2fc-7) 



(2) 



where n(r) = njsM for the constant density case and n(r) 
for the wind environment with A = J^y g/cm. 



, 2 2 (^''~^^) 

e — 4nmpC 7^1^(12-3(1)^ 



and 

7 = 7/X" 



(3) 



(4) 



To solve for the above quantities, we need to find 7/, the 
Lorentz factor of the material directly behind the shock front 
that defines the blastwave outermost edge. It is given by 
(e.g.. Chevalier & Li, 2000) 



If = 



{17-4k)E 



l6'Kmpn{R)R-^c^ 



(5) 



where E is the energy of the blastwave. 

Let us consider now the point P in Fig. [1] It is charac- 
terized by a distance r from the center of the explosion and 
an angle 8 with respect to the line of sight. Velocity at this 
point is in the r direction. The point P begins to emit radia- 
tion when it is crossed by the forward shock at the observed 
time: 



Tabs ^ ~ ^ ^ 



cos{e) + 



4(4 - k)'yj 



(6) 



The radiation observed at time Tots is obtained therefore by 
integration over the locations (r, 9) that satisfy: 



^11 



cos(9) 



1 



4(4 - kWj 



< Tot 



(7) 



At any location that satisfies Eq. [T] the ratio R/r that 
is needed to compute the self-similar variable x is computed 
through the equation: 



Tob 



R 

c 



4(4 ~ khj 



■ cos{e) 



(8) 



where R is the radius of the shock at the laboratory time at 
which the photon was emitted at radius r to be observed at 
infinity at time Tots- 

The energy density of the magnetic field and of non- 
relativistic particles is assumed to be a constant fraction 
of the local internal energy, as derived from the BM blast 
wave solution. Such a prescription is analogous to adiabatic 
evolution of the population of relativistic electrons with the 
density of particles (Beloborodov 2005). We use the full form 
for Weibel fields described in detail in Medvedev (2006) and 
Workman et al. (2008) to describe the geometry and cor- 
relation function of the magnetic field. We assume that the 
properties of the magnetic field stay constant throughout the 
afterglow and throughout the blast wave. We parametrize 
the magnitude of the magnetic field according to the usual 
prescription 



{Bf 
8-K 



— ese 



(9) 



where en represents the equipartional fraction of the mag- 
netic field with the energy density. 



2.2 The Observed Jitter Spectrum 

To determine Plat{ij), the total power emitted per unit fre- 
quency by an ensemble of electrons radiating in a unit vol- 
ume, we perform the following integral 
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PLtH = 



n {'-f)p' {uj)d'y. 



(10) 



The number density n and the minimum Lorentz factor 
{■ymin) in the integral are derived by basic assumptions 
about the total number and energy of electrons, the inter- 
ested reader is referred to Workman et al. (2008). Once the 
above quantity has been computed we assume isotropic radi- 
ation and, following Rybicki and Lightman (1979), we define 
the emission coefficient, as 



Ju. 



47r 



(11) 



Fofiowing the process outlined in Medvedev et al. (2007), 
Workman et al. (2008); Granot, Piran, and Sari (1999b), 
and Rybicki and Lightman (1979) we define the absorption 
coefficient as 



7 



-(p+i) p' 



P'{u))d'y. 



(12) 



Once we have found the necessary coefficients we use the 
fact that UloUJ and ^ are both Lorentz invariants with 
uj = aj'76(l -I- -^cosS'); where 75 and v refer to the bulk 
Lorentz factor and velocity of the afterglow and cos6'' is 
defined by the angle between a photon's path and the co- 
moving volume's direction. Using this property we can re- 
write the differential equation for specific intensity in our 
frame, given by 

= jui - oi^Iu) (13) 

directly in terms of comoving quantities. We assume that 
there is no scattering or frequency redistribution (see Fleish- 
man 2006ab). This allows us to calculate the spectrum by 
numerically solving for a series of one dimensional, line of 
sight, integrals. We define our integration region as the three 
dimensional volume of constant arrival time, which is given 
by equation [S] As the jitter radiation depends on the angle 
of emission in the co-moving frame, we solve for 6' using 6 
and /3 by inverting the equation 

cos{e') + p 



cos{e) = 



1 -I- pcos{e') 



(14) 



Once we have found the comoving angle at each point, we 
use the the Blandford-McKee solution to find the values of 
7, n', e', and the magnetic field at each point and solve for 
ji and a'^. 

With these functions known throughout the volume of 
constant arrival time, we solve the equation for observed 
specific intensity, at a given frequency, along lines of sight 
that span the vertical extent of the surface, giving us I^jiy) 
where y is the perpendicular distance from the center of the 
source (see Fig. [1} . 

Cooling is added analytically by changing the spectral 
slope to uj^~'^ above a frequency of 

ll6 



8.5 X 10'- ^-1/2/ 



3/2^-1 ^-1/2 
ISJ\I days 



27V '0.001 
(from Sari et al. 1998) for the ISM case and 

5.7 X 10^^ 



27r 



_^-i/2._e£_s- 
^0.001^ 



3/2 --^2,1/2 
^* ''days 



(15) 



(16) 



(from Chavalier & Li 2000) for the wind case, where A* — 
■ixio^^ cm-i • However, for the parameters and times consid- 



ered in this paper, cooling only becomes important above 
32 keV for the ISM case and 24 keV for the wind case. 

Once we have solved for the specific intensity in the 
observer frame, the observed flux at the detector is given by 



= 2tt- 



d? 



(17) 



where d is the distance to the afterglow and R_L,max is the 
largest perpendicular extent of the afterglow (as in Granot 
et al. 1999a). 



3 RESULTS 

In order to compare synchrotron and jitter radiation we cal- 
culate the spectrum and time evolution for both radiation 
mechanisms for the case of a constant (ISM) and a wind ex- 
ternal medium, for a total of 4 models. Analytical lightcurves 
for these models are derived in Medvedev et al. (2007) . In all 
four cases the total isotropic equivalent energy in the blast 
wave was E = 10^'' erg, the magnetic field equipartition 
fraction was ~ 10~^, and the fraction of internal energy 
in electrons was — 0.01. For the ISM case the external 
density was n — 1 while for the wind medium we assumed 
M = lO"^M0/yr and = 1000 km/s giving A, = 0.1. 

A spectrum for each of these 4 models ai t — 500s is 
plotted in Fig. [51 At this time, the two synchrotron models 
have a peak frequency of about 0.3 eV. At lower frequencies 
the spectrum increases as oj^^^ and at higher frequencies 
it decreases as oj"'''"^'/^. The self-absorption break is ob- 
served at ~ lO"'^ eV. The two jitter radiation models have 
similar behavior to the synchrotron models above the peak 
frequency and in the optically thick regime. However, the 
jitter models have a nearly flat spectrum below the peak 
frequency rather than a ly^''^ branch. Even though jitter ra- 
diation can produce a spectral slope as steep as lo^ at small 
angles, the integrated emission is dominated by material at 
large comoving angles (see Fig. [3} which produce a flat spec- 
trum for jitter radiation (see Medvedev 2006; Medvedev et 
al. 2007; Workman et. al. 2008). A good quality spectrum 
over a wide frequency range can distinguish between the syn- 
chrotron and jitter radiation mechanisms if the v^^^ or fiat 
portion of the spectrum is captured. However, the spectrum 
alone cannot separate ISM and wind external media. 

At very low frequency (see Fig.|3)) all models are initially 
optically thick. This results in an increase of luminosity as 
t^^'^ for the two ISM models and an increase as for the 
two wind models. The radiation mechanism does not signif- 
icantly afi'ect the results in the optically thick regime. 

At higher frequency all four models are optically thin. A 
comparison of the light curves for the four models at 1 eV is 
shown in Fig. [5] For the ISM case, the synchrotron spectrum 
initially increases as t^^'^ and then transitions to a steep 
decay (f'^^fP"!)/*) as the peak frequency of the spectrum 
decreases below 1 eV. This light curve is similar to the low 
frequency synchrotron-ISM and jitter-ISM models is Fig. [l] 
With a wind external medium and synchrotron radiation, 
the light curve is initially flat before transitioning to a steep 
decay {t~^^^~^^^'^). If the jitter radiation mechanism is used, 
the slope of the light curves are initially decreased by a factor 
of t~^''^, giving a flat initial slope for the jitter-ISM model 
and a decrease of t~^^^ for the jitter-wind case. The result 
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is that in this regime the jitter-ISM and synchrotron-wind 
models are not distinguishable from each other from the 
light curve alone, as they both begin with a flat intensity. 

For the jitter-wind case, the luminosity initially de- 
creases as t~^^^. This gives the light curve the appearance 
of a broken powerlaw decay. The break occurs at about lO^s 
at 1 eV for this example, but the time of the break depends 
on frequency as ttreak oc uj~^^^. Fig. |6] shows light curves 
from the jitter-wind model for energies between 10~^ eV 
and 1 keV. The jitter wind model naturally gives rise to a 
power law decay with a chromatic break between slopes of 
t"^/^ and t-(3p-i)/4 At lower frequencies in Fig.[6l the ini- 
tial slope is increasing as the ejecta are initially optically 
thick. Observations of the jitter-wind model would show a 
chromatic break over a range of frequencies or a steep decay 
at high frequency with a shallow decay and then a break at 
lower frequencies, depending on the time range and frequen- 
cies of the observations. 

Fig. [7] shows light curves from the jitter- wind model in 
the x-ray (1 keV), visible (575 nm), and radio (24.2 GHz). 
In this case, the x-ray observations show a single powerlaw 
decay, the visible shows a broken powerlaw, and the radio 
shows a powerlaw increase. 



4 DISCUSSION AND CONCLUSIONS 

We have calculated the spectra and light curves produced 
by jitter radiation from a relativistic fireball expanding into 
either an ISM (constant density) or wind (p (x r~^) exter- 
nal medium including self absorption and electron cooling, 
and then compared these results to results for synchrotron 
radiation in identical cases. We improved upon our previ- 
ous computations (Medvedev et al. 2007) by performing the 
integration over the blastwave and taking into account the 
light travel time of the photons. 

The primary difference between jitter and synchrotron 
radiation is that jitter radiation produces a flat spectral 
slope below the peak frequency, rather than the uj^^^ slope 
seen in synchrotron radiation. Jitter and synchrotron radia- 
tion are not distinguishable in the optically thick regime or 
above the peak frequency, only in between. The flat spec- 
tral slope of jitter radiation has the effect of decreasing the 
slope of the light curve by a factor of t~^^'^ compared to 
synchrotron radiation in the relevant frequency and tempo- 
ral intervals. This means that the slope of the light curve, 
while the observed frequency is below the peak frequency, 
decreases from t^^^ for synchrotron to for jitter in the ISM 
case and decreases from to t~^^'^ for the wind case. This 
indicates that the jitter-ISM and synchrotron-wind cases 
would not be distinguishable from a single light curve as 
they both have a slope. 

In the last four years, Swtft observations, combined 
with ground-based optical and radio telescopes have put the 
standard synchrotron external shock model under scrutiny. 
The new observations have shown temporally flat X-ray 
light curves (Nousek et al. 2006), chromatic breaks (e.g., 
Huang et al. 2007), and uncorrelated multi-band behavior 
(e.g., Blustin et al. 2006; Dai et al. 2007). The light curves 
and spectra of jitter and synchrotron radiation are differ- 
ent. However, the differences appear to be too subtle to ex- 
plain the atypical afterglow observations summarized above. 



Even though jitter radiation can produce flat lightcurves, it 
does so in all wavebands and not only in the X-rays, and 
cannot explain the steep decay observed, for example, in 
GRB 070110 (Troja et al. 2007). Chromatic breaks can be 
produced by jitter radiation as well as by synchrotron, the 
only difference being the amount of change in the pre- and 
post-break decays. Pending a more quantitative compari- 
son of the jitter light curves and spectra with the data, it 
seems that jitter is a viable radiation mechanism to explain 
the behavior of GRB afterglows. It shares, however, all the 
problems of synchrotron in the explanation of some of the 
most recent data. Some possible explanation and a review of 
the difficulties can be found, for example, in Ghisellini et al. 
(2007); Liang, Zhang & Zhang (2007); Uhm & Beloborodov 
(2007); Panaitescu (2008). 
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The line A-B represents a line of sight along which the radiative transfer equation is solved to find Ito{y)- 
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Afterglow Spectrum at 500s 
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Figure 2. Comparison of afterglow spectra at t = 500 seconds. Spectra plotted are for radiation from a shock expanding into a wind or 
ISM-type medium produced by jitter and synchrotron radiation. Both jitter spectra are flat between the peaJt frequency and optically 
thick cutoff, while the synchrotron spectra increase as ui^/^ before reaching a broad peak. The wind cases are optically thick at lower 
frequencies. 
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Surface Brightness 
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Figure 3. Comparison of surface brightness at 1 eV at t = 1000 seconds for the same 4 models of Fig. [2] In all cases, the emission is 
concentrated towards the outer edge. Normalized transverse distance is defined as y/ymax with y defined as in Fig. [T] 
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Time Evolution at 10"^ eV 
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Figure 4. Simulated light curves at hu = 10"^ eV (A = 1.24 mm) for the 4 models of Fig. [2] All 4 models are optically thick at this 
frequency. The ISM models decrease in luminosity as t^/^, while the wind models increase as . 
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Time Evolution at 1 eV 
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Figure 5. Simulated liglit curve at hu = I eV (A = 1240 nm) for the 4 models of Fig. [2] Tlie synchrotron-ISM model luminosity increases 
at early times as t^^^. The synchrotron- wind and jitter-ISM models are both initially flat before transitioning to a power-law decay as the 
peak frequency decreases below hu = 1 eV. The jitter- wind model initially decreases as t~^^^ before transitioning to a steeper power-law 
decay. 
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Time Evolution Jitter-Wind nnodel fronn 1 eV to 1 keV 
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Figure 6. Comparison of light curves for jitter-wind model at different energies. The transition time between the luminosity scaling as 
t~^/^ and t-C^P"!)/* is dependent on frequency and gives the appearance of a chromatic break in a power-law decay light curve. The 
time of this brealc ranges from about t = 10s to t = lO^s for the frequencies shown in the figure. 




Figure 7. Comparison of light curves for jitter-wind model in the X-ray {hf = 1 keV), visible (A = 575 nm) and radio {i/ = 24.2 GHz) 
bands, observed from t = 60s to t = lO^s. The x-ray light curve is a single powerlaw decay, the visible is a broken powerlaw and the 
radio, in the optically thick regime, is a powerlaw increase. 



